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NOTATION 


FORTRAN 

Symbol Name Description 

Number  of  stations. 


N 

NO  (I) 

x.  DISTL  (I),  A 

d DIA  (I),  B 

a 

(a)  STIFX(I),  C 

(b)  I El  or  ST 

(a)  dc  WTS(I),  D 

(b)  T Ap 

W CONWT(I),  E 

a IM  (I),  I 

b JM  (I),  J 

IBRG(I),  K 
IS  (I) 


Station  number  of  I1*1  station. 

Length  to  extreme  end  from  I1*1 
station  in  inches. 

Outer  diameter  of  shaft  at  I1*1 
station  in  inches. 

(a)  Inner  (if  hollow)  or  second  (if 
sleeve)  diameter  in  inches. 

(b)  Stiffness 

(a)  Third  diameter  in  inches 

(b)  Weight  in  lb. /in.  per  unit 
length, 

(a)  Concentrated  Weights  in  lbs 
(negative  in  sign). 

(b)  and  Bearing  reactions  (positive 
in  sign). 

Outer  material  at  I*h  station 
(only  1~5  possible). 

Inner  material  at  I1*1  station 
(only  1-5  possible). 

Code  indicating  presence  (1)  or 
absence  (0)  of  bearing  at  station  I. 

Code  indicating  sand  (1)  or  no 
sand  (0)  in  bore. 


FORTRAN 


Symbol 

Name 

Description 

E 

YMAT,  H 

Moduli  of  elasticity  in  lbs// sq.  in. 
(Young's  modulus). 

$ 

DENMAT,  P 

Density  of  shaft  materials  in  lbs.  / 
cu.  in. 

li  * xi+rxi 

DIST  T(I) 

(a)  Section  length  in  inches. 

(b)  Stiffness  (10  10/EI) 

ML  - / o SdL 

BMTAB 

-3 

Bending  moment  in  10  inch-lbs. 

e = i / l MdL 
El 

THETA 

Slope  in  Radians. 

6L  = r 0 eL  dL 

DEF 

Deflection  in  inches. 

Of 

BNFLU 

Influence  numbers  for  bearings 
(lb.  per  . 001- inch  rise)  and  de- 
flection and  slope  (for  1-inch 
rise). 

DIBRG 

Sec f ion  length  from  support  point 
to  extreme  end  table. 

IB  TAB 

Bearing  index  table.  v 

n 

IBRS 

Number  of  bearings. 

SL  = J'oWdL 

SFTAB 

_ 3 

Shear  force  in  10  lbs. 

IERROR  Error  code  in  incorrect  data 

preparation. 


Symbol 


FORTRAN 

Name 


Description 


TITL2 

IXM 

SWAT 

NS 

IYM 

Case  title  array  for  MGE2. 
Sentinel 

) Unused  options. 

R 

CANOT 

CONOT 

Bearing  Reaction  in  lbs. 

WBAR 

Total  Shaft  weight  including  con- 
centrated loads  in  lbs. 

XBAR 

Distance  of  C e nte r-of-g r avity 
from  station  no.  1 in  inches. 

A 

a 

AREA  1 

Shaft  cross  section  area,  second 
material. 

*b 

AREA  2 

Shaft  cross  section  area,  second 
material. 

A 

c 

AREA  3 

Shaft  cross  section  area,  sand  in 
bore. 

'a 

AMI1 

Moment  of  Inertia  for  first 
material. 

h, 

AMI2 

Moment  if  Inertia  for  second 
material. 

XBRG 

Bearing  relocation  off  straight 
line. 

TITL5 

Case  title  array  for  MGE5. 

Some  storage  areas  have  uses  a and  b in  different  sections  of 
the  program.  The  several  arrays  which  have  special  names  in  READIN 
subroutine  are  listed  as  second  entry  in  name  column. 
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ABSTRACT 

MGE2  is  a FORTRAN  IV  program  which  computes  characteristics 
of  the  main  propulsion  shaft  of  a ship,  assuming  that  all  bearine 
centers  are  on  a straight  line.  Output  from  this  program  consists  of 
the  straight-line  bearing  reactions  and  the  influence  numbers  for  each 
support  point.  These  influence  numbers  indicate  the  changes  in  supporting 
reactions  produced  by  the  raising  or  lowering  of  a particular  bearing. 

The  program  also  calculates  the  shear,  bending  moment,  slope,  and 
deflection  values,  and  the  weight  and  stiffness  factors  for  stations  desig- 
nated along  the  shaft. 

MGE5  is  a program  similar  to  MGE2.  It  is  a design  tool  for 
determining  the  effects  of  bearing  centers  at  any  given  elevation  relative 
to  the  straight  line.  MGE5  is  often  used  after  MGE2  to  determine  the 
reactions  of  realigned  bearings.  Output  from  MGE5  includes  bearing 


reactions  at  each  support  point  and  shear,  bending  moment,  and  slope 
and  deflections  values  at  designated  stations. 

\ 

ADMINISTRATIVE  INFORMATION 
The  work  reported  herein  was  authorized  under  job  number 


1 

' I 

INTRODUCTION 

MGE2  gives  the  characteristics  of  a shaft  with  all  hearings, 
or  support  points,  on  a straight  line.  Output  includes  shear,  bending 
moment,  slope,  and  deflection  as  well  as  weight  and  stiffness  factors 
at  designated  stations  along  the  shaft.  A table  of  influence  numbers  is 
listed  for  each  bearing.  These  influence  numbers  indicate  the  changes 
in  supporting  reactions  caused  by  the  raising  or  lowering  of  a particular 
bearing. 

MGE5  calculates  shear,  moment,  slope,  and  deflection  at  various 
stations  and  bearing  reactions,  with  bearing  centers  set  in  any  given 
elevation  relative  to  a straight  line. 

The  programs  were  coded  in  FORTRAN  IV  with  the  exception  of 
MAP  double  precision  routines  in  the  7090  matrix  inversion  portion. 

Decks  or  further  information  are  available  from  Code  830,  Naval 
Ship  Research  and  Development  Center  (NSRDC),  Washington,  D.  C. 

20007. 

These  programs  have  also  been  run  on  LARC  and  7090  in 

1 .1 

FORTRAN  II. 


METHOD  OF  CALCULATION 


(Taken  from  the  "Calculation  of  Ship  Propulsion  Shafting  Bearing 
Reactions  on  an  IBM  650  Computer"  written  by  Edward  T.  Antkowiak) 

"Although  the  texts  on  indeterminate  structures  illustrate  various 
devices  for  handling  continuous  beams  on  multiple  supports,  such  as 
Hardy  Cross  or  Relaxation  Method,  Three  Moment  Equations,  etc. , 
the  large  memory  capacity  and  speed  of  a digital  computer  permits  a 
more  'classical'  approach.  This  approach  is  based  on  a direct  evalu- 
ation of  the  relationship  of  the  beam  lead  (w)  to  the  deflection  (6). 


"Therefore,  successive  integrations  of  the  load- stiffness  diagram 
will  yield  shear,  moment,  slope  and  deflection.  However,  in  order  to 


obtain  zero  deflection  at  each  support,  the  correct  reactions  would  have 
to  be  included  in  the  load  diagram  prior  to  the  integration.  In  the 
shafting  problem,  the  reactions  at  each  support  are  the  unknowns. 

"If  the  reactions  are  not  included  in  the  load  diagram  and  it  is 
assumed  that  the  beam  is  held  firmly  at  one  end,  the  successive 


integrations  will  result  in  a deflection  diagram  similar  to  that  shown 


in  Figure  1 (a).  Applying  estimated  reactions  whose  magnitudes  are 
proportional  to  the  deflection  will  yield  a diagram  with  a smaller  de- 
flection, as  shown  in  Figure  1 (b).  Further  approximations  and 
sucessive  integrations  will  converge  to  a zero  deflection  at  the 


supports,  as  shown  in  Figure  1 (c).  At  this  point,  the  correct 


reaction  has  been  applied 


Figure  1 - Convergence  of  Shaft  Deflection  Curve  to  a 
Straight  Line  by  Adding  Support  Reactions 


In  order  to  apply  this  principle  to  practical  use,  some  method  of 


obtaining  reaction  estimates  had  to  be  devised.  The  method  used  is  to 


determine  the  amount  of  reaction  change  produced  by  raising  or  lowering 


a bearing  support  a unit  amount.  This  is  also  one  of  the  final  answers  to 
our  problem.  These  influence  numbers  are  obtained  by  assuming  the 
shaft  has  no  weight  and  that  a unit  reaction  is  applied  to  one  bearing. 
Successive  integrations  will  then  yield  resulting  deflections  at  all  other 
bearing  locations.  By  applying  this  unit  reaction  to  each  support  in  turn, 
a family  of  deflections  will  be  obtained  which  can  be  substituted  into  a set 


Vi* 


of  equations  (two  more  equations  than  there  are  number  of  support  points), 
expressing  the  relationship  of  reactions  and  deflections.  Solution  of 
these  equations  will  yield  the  'influence  numbers'  necessary  to  the  final 
solution  of  straight  line  reactions. 

"Multiplying  the  influence  numbers  by  the  deflections  obtained  from 
integration  of  the  unsupported,  weighted  (real)  shaft,  and  summing  the 
reaction  changes  for  each  bearing,  the  final  straight  line  reactions  will 
result.  Applying  this  total  reaction  for  each  bearing  to  the  load  diagram 
and  repeating  the  integration  process,  the  deflection  at  each  support  will 
be  zero,  thus  proving  that  the  correct  reaction  has  been  applied.  If  the 
deflection  is  not  zero,  the  entire  process  can  be  repeated  until  the 
bearing  deflection  becomes  trivial.  " 
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I HE  OR 


Figure  2 shows  a sample  section  of  shaft  with  nomenclature  as  in 
the  Notation.  The  programs  compute  moment  of  inertia  and  weight  per 
unit  length  as  follows: 


Outer  area 


Inner  area 


Innermost  area 


Inertia  of  outer  part 


Inertia  of  inner  part 


Weight  per  unit  length 


Figure  2 - Section  of  Shaft  with  Computational  Inputs  Labelled 


where  the  first  term  (Material  a ) must  be  present  in  all  cases;  the 
second  term  (b)  will  be  present  for  two  materials;  the  third  (c)  is 
present  only  if  the  bore  is  filled  with  sand;  and  the  fourth  (d),  which 
represents  the  upthrust  of  water,  is  present  only  if  the  shaft  is  wholly 
or  partly  immersed  in  water. 

Consider  a shaft  of  no  weight  but  with  the  same  stiffness  as  the 
real  weighted  shaft.  If  such  a weightless  shaft  were  on  "straight  line" 
bearings,  there  would  be  no  reaction;  however,  if  one  of  the  bearings 
should  be  raised  by  1 mil,  forces  would  occur  at  each  bearing. 

These  forces  are  the  "influence  numbers"  which  form  an  n x n matrix 
where  n is  the  number  of  bearings.  The  element  in  the  ith  column  and 
jthrow  gives  the  force  on  Bearing  i for  a deflection  of  1 mil  at  Bearing  j. 

Consider  the  weightless  shaft  acting  as  a cantilever  as  shown  in 
(a)  of  Figure  3.  When  a unit  force  is  applied  at  B,  Section  AB  will  be 


straight,  but  inclined  at  an  angle  to  the  horizontal  (b).  There  is  an 


advantage  to  assuming  that  Section  AB  remains  horizontal  with  the 
wall  moved  (c),  using  AB  as  the  reference  axis. 


Now  consider  the  shear  force,  bending  moment,  slope,  and 
deflection  diagrams.  According  to  the  theory  of  beam  deflection,  these 
diagrams  are  related  to  one  another  by  integration  (see  Figure  4). 


Figure  4 - Successive  Integrations  of  a Load  Stiffness  Diagram 


The  result  of  a series  of  such  integrations  is  summarized  in  tabular 


form  by  Table  1.  Thus  d„.  is  the  deflection  at  Station  4 for  a force 


of  unity  applied  at  Station  2 


TABLE  I 


If  the  beam  Section  AB  is  not  considered  to  be  horizontal  or  at 


zero  slope,  these  deflections  must  be  corrected  as  follows: 

At  A:  add  Aq 

At  B:  add  A + fi  x. 

o o 1 

At  C:  add  A + 9 x„  (see  Figure  5). 

o o 2 


4 original  bed»r  I 


Figure  5 - Beam  (Shaft)  Displaced  from 
Horizontal  Reference  Line 


When  unit  reactions  Rj,  R£,  ...  Rfi  occur  at  support  points  1,  2,  ...  n 
along  the  weightless  shaft,  the  following  equations  apply: 


6 = A + 9 x.. 

1 o o 1 


VV  flox2  * d12  R1 


6 = A + 0 + x + d1  R*  + d„  + . . . d, 
n o o n In  1 2n 


R/ 


(n-l)n  (n-1) 

Two  additional  equations  are  required  to  solve  this  group;  they  may  be 
found  by  considering  the  beam  end  C in  Figure  5 as  unfixed  when  it  is 
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supported  by  reactions  along  its  length.  When  acting  as  a multiply 
supported  beam,  this  free  end  has  no  shear  or  moment,  such  that 


0 — R,  + R_  + . . . + R 
12  n 


0 = x.  R,  + x„  R„  + . . . x R 
112  2 n n 


All  these  equations  can  therefore  be  put  into  a n+2  equation  square 
matrix  form 


.>  v_ 

1 XjO  0 

0 ... 

i 

o 1 

62 

lX2d12  ° 

0 ... 

. 

% , 

63 

• 

1 = 

lX3  d13  d23 
• • • • 

0 ... 

* • • • 

X 

Ri 

6 

n 

. . . 

O • • • 

IX  d1  d„ 
n In  2n 

• • • • 

• • • • 

d3n 

0 

0 0 1 1 

1 ... 

• 

0 

-\  r 

o o xi  x2 

\ 

R 

s n. 

representing  the  deflections  for  a given  reaction.  Inversion  of  this 
matrix  gives  the  reactions  for  a given  deflection,  the  reaction  influence 
numbers.  The  columns  which  give  deflection  and  slope  condition  now 
give  slope  and  deflection  influence  numbers  (i.  e. , the  end  slope  and 
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the  deflection  for  a given  deflection  at  an  intermediate  point).  Thus  a _ 

OCt 

represents  the  reaction  at  Support  3 for  a theoretical  unit  deflection 
at  Bearing  2. 


Assuming  the  initial  slope  and  deflection  are  zero  and  that  reactions 
at  support  points  are  zero  (shaft  cantilevered  at  one  end),  the  end  shear 
(Vj)  and  end  moment  (M^)  of  the  weighted  shaft  are  calculated.  These 
results  are  used  to  modify  the  first  (R^)  and  last  (R^)  reactions,  so  that 
Vj  and  then  approach  zero. 

Consider  the  actual  shaft  with  only  two  external  bearings,  held 
at  the  end  to  ensure  no  initial  slope  and  deflection  (see  Figure  6),  with 
total  shaft  weight  W,  total  moment  from  the  other  end  M,  and  total 


The  above  discussion  presents  the  solution  to  the  two-bearing 
problem.  Then,  by  using  the  influence  numbers,  the  effect  on  the 
reactions  when  all  the  bearings  are  moved  to  the  correct  position  can 
be  found;  Aq  and  0Q  are  modified  at  this  time. 

The  actual  answers  from  the  inversion  will  not  be  exact,  but 
repeated  entry  into  the  deflection  integration  routines  will  give  corrected 
values  of  the  deflection  for  all  bearings.  Several  of  these  weighted 
integrations  may  be  required  to  obtain  results  within  the  required 


margin 


MGE2  - STRAIGHT-LINE  BEARINGS 


PROGRAM  DESCRIPTION 

To  calculate  the  characteristics  of  a shaft  with  all  supports 
on  a straight  line  as  outlined  in  "Method  of  Calculation,"  the  shafting 
system  must  be  divided  into  a number  of  stations,  subject  to  the 
following  restrictions: 

1 Total  number  of  stations  must  not  exceed  400. 

2 Total  number  of  supports  must  not  exceed  22 

3.  Each  change  in  cross-section  diameter,  each  bearing  location, 
each  concentrated  load,  and  the  beginning  of  each  immersion  requires 
a new  station 

4 Although  the  length  of  each  station  is  arbitrary,  there  should 
be  enough  of  these  stations  to  provide  sufficient  abscissa  points  for 
plotting  shear,  moment,  slope,  and  deflection  diagrams,  if  desired. 

5 All  stations  (except  the  last)  must  have  nonzero  diameter  and 
material  The  last  data  card  is  not  considered  a section;  it  serves  to 
specify  the  total  length  of  the  shafting  system. 

6 Bearing  may  not  occur  at  the  first  or  at  the  last  station. 

7 Concentrated  loads  should  not  appear  at  bearing  locations 
since  calculated  bearing- reaction  values  will  overwrite  the  concentrated 
loads  Add  an  additional  shaft  section  with  a 0.  100- inch,  or  even  a 
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0 0001 -inch,  separation 


I 


The  program  is  written  in  a general  manner  and  provides  for  the 
following  variables: 

1 Material  of  shaft 

2 Shaft  in  air  or  water 

3 Shaft  solid  or  hollow 

4 Sand  in  bore  of  hollow  shaft 

The  specific  moduli  and  densities  in  the  program  are: 

Code  Material  Density  Modulus 
(lb/cu  in. ) (lb/sq  in. ) 


1 

Steel 

-.28355 

30  x 106 

2 

Bronze 

-.  31399 

15  x 10* 

6 

Sand 

-.0636 

7 

Salt  water 

-.03705 

If  materials  used  in  a specific  shaft  have  constants  quite  different 
from  the  above,  additional  values  may  be  added 

Slight  changes  in  diameter  may  be  ignored  when  the  effect  on 
weight  and  stiffness  of  the  shaft  is  negligible  The  average  diameter 
of  the  propeller  taper  is  used. 

Loads  acting  downward  will  be  negative.  In  most  cases,  reactions 
will  be  positive.  However,  in  some  instances  where  several  bearings 
are  close  together,  one  bearing  may  indicate  a negative  reactioafor 
the  straight-line  condition 

The  weight  of  the  propeller  is  applied  as  a concentrated  load  at 
its  center  of  gravity,  if  known;  otherwise,  the  weight  is  considered  to 
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act  at  the  middle  of  the  tapered  section,  The  dry  weight  of  an  immersed 
propeller  may  be  reduced  to  take  into  account  the  buoyant  effect  of 
water  by  using  the  relationship 

D -D 

W (water)  - W (air)  * — O — — 

P 

where  D is  the  density  of  the  propeller  material  and  D is  the  density 
p w 

of  the  water. 

Bearing  reactions  are  assumed  to  act  at  points  in  the  middle  of  the 
fore  and  aft  ends  of  the  bearings.,  However,  in  the  case  of  a relatively 
long  bearing  such  as  a stern  tube,  it  may  be  assumed  that  the  load  is 
concentrated  one  shaft  diameter  forward  of  the  aft  end  of  the  bearing. 

For  a more  detailed  study,  support  points  at  several  positions  along  a 
bearing  may  be  assumed. 

The  program  gives  the  magnitude  of  all  bearing  reactions,  assuming 
all  bearings  on  a straight  line,  It  prints  slope  and  deflection  influence 
numbers,  which  represent  the  amount  of  change  in  all  bearing  loads 
resulting  from  raising  a particular  bearing  0,  001  in.  from  the 
straight  line.  Values  of  shear,  moment,  slope,  delfection,  weight 
(in  pounds  per  inches),  and  stiffness  (in  terms  of  -gj)  are  given  for 
each  station.  The  total  weight  of  the  shaft  and  the  location  of  the  center 
of  gravity  are  also  printed. 


I 
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The  criteria  for  determining  zero  deflection  at  a bearing  and  the 
absence  of  shear  force  or  moment  at  the  free  end  are  listed  below: 

1.  Deflection  at  a bearing  no  greater  than  0.  001  in. 

2.  Shear  force  at  free  end  less  than  10  lb. 

3.  Bending  moment  at  free  end  less  than  100  in.-lb. 

4.  A maximum  of  four  iterations  on  the  deflection  check. 

5.  A maximum  of  four  iterations  on  the  shear  and  moment  checks. 


The  following  error  printouts  occur  after  unrecoverable  errors 
have  been  detected.  The  program  then  transfers  to  next  data  case. 
The  printout  TOO  MANY  ITERATIONS  DEFL  XXXX  MOMENT 


are  not  input  in  ascending  order  of  distance  from  Station  1.  Case  is 
terminated. 

The  printout  SINGULAR  INFLUENCE  NUMBERS  occurs  if  matrix 
is  singular.  Since  it  is  indicative  of  incompatible  data,  case  is 
terminated 


XXXX  occurs  when  four  iterations  through  DEFLN  and  BEAM  have 
failed  to  meet  terminal  deflection,  bending  moment,  and  shear  force 
criteria.  The  current  results  are  then  accepted. 

The  printout  STATION  DISTANCE  ERROR  occurs  when  stations 


OPERATING  INSTRUCTIONS 


1.  The  program  run  with  7090  IBSYS  Version  13  includes  the 
following  cards: 

$JOB 

^EXECUTE  IBJOB 

$IBJOB  FIOCS 

decks  MGE2 
RED2 
BEAM 
DEFL 
INOU 
INFN 
AM  MAT2 
AM  DPAF 

$DATA 

MGE2  heading  card  for  input  case  1 
One  card  per  station  for  input  case  1 
One  card  for  each  constant  change  (omit  if  no  material 
constant  changes) 

Sentinel  card  with  punch  in  column  2 
Additional  cases  of  data,  if  called  for  by  sentinel 


Final  sentinel  has  9-  or  7-punch  in  column  2 
l END  OF  FILE 

2.  Output  of  five  or  more  pages  per  case  is  written  on  SYS0U1  for 
program  control  printing. 

At  NSRDC,when  compiled  with  plot  option,  control  cards. 

fATTACH  B9 

$ AS  SYSLB4 

are  inserted  before  $ EXECUTE  card,  tuid  plots  for  SC  4020  may  be 
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prepared  on  SYSLB4.  Deck  GRAP  must  also  be  included  before 
$DATA  card. 

3.  Total  time  for  IBSYS  compile,  load,  and  execute  of  six  varied  cases 
was  1. 8 min.  Since  the  load  time  for  retrieval  of  binary  program 
deck  from  master  instruction  reel  with  $IEDIT  is  1 min,  it  is  obviously 
advantageous  to  run  multiple  cases,, 

4.  No  sense  switches,  sense  lights,  or  scratch  files  are  used  by  the 
program.  Memory  requirement  (including  I/O  subroutines)  19,400. 

For  smaller  computer,  dimensions  could  be  revised. 

5.  The  program  run  on  LARC  includes: 

Data  cards  for  Program  MGE2,  followed  by  six  cards  with  END  OF 
TAPE  punched  in  columns  1-12,  are  converted  to  tape  and  ready  for 
mounting  on  Unit  11. 

Object  code  tape  containing  instructions  for  routine 

MGE2 

READ2 

BEAM 

DEFL 

INOU 

INFN 

AM  MAT2L 
is  mounted  on  Unit  13. 

Output  is  written  on  Tape  20  for  l/l  LSC  loop  printing. 
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SC  4020  PLOTS 


An  option  for  SC  4020  graphic  output  may  be  used  at  NSRDC  on 
7090  or  LARC.  Card  MGE2  0164  (or  MGE5  0119)  is  changed  from 
981  CONTINUE  to  981  CALL  GRAPH 

The  graph  routine  is  listed  in  Appendix  A. 

The  two  sets  of  plots  made  for  each  case  are  moment  and  deflection 
versus  station  distance.  If  the  total  number  of  stations  is  less  than  99.  one 
frame  is  produced  for  each  graph.  For  larger  numbers  of  stations,  a 
fixed  station  distance  scale  of  450  ft/frame  results  in  2 to  21  frames 
for  each  graph. 

The  ordinate  versus  abscissa  is  plotted  as  a continuous  line  with 
bearing  stations  further  identified  by  X.  See  Appendix  A for  sample 
output. 
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SHAFT  ARRANGEMENT 


Figure  7 shows  a typical  line  shaft  dimensioned  for  preparation  of 
input  data.  The  sketch  should  be  made  to  some  appropriate  scale,  using 
information  from  arrangement,  assembly,  and  detailed  drawings  of  low 
speed  gear  and  shaft.  The  relationship  of  the  components  can  then  be 
readily  observed. 

1.  The  total  number  of  stations  must  not  exceed  400  and  the  total 
number  of  support  points  must  not  exceed  22. 

2.  Station  1 is  the  forward  end  of  the  shaft  and  all  other  stations 
are  located  in  terms  of  their  distance  from  Station  1.  A station  must  be 
initiated  at  each  support  point,  at  each  concentrated  load,  and  at  each 
significant  change  in  cross-section  diameter.  Concentrated  loads 
should  not  appear  at  bearing  locations;  an  additional  station  with  0.  001-  or 
0.  0001-in.  separation  may  be  used  (see  note  at  end  of  paragraph).  To 
reduce  the  number  of  stations,  slight  changes  in  diameters  may  be 
ignored  or  combined  when  their  effect  on  the  weight  and  stiffness  of  the 
shaft  is  negligible.  Thus,  in  Figure  7,  since  the  journal  sections  were 
only  1/4  in.  wider  than  the  shaft,  they  were  considered  to  be  the  same 
diameter  as  the  shaft.  Also  the  average  diameter  of  the  tapered  propeller 
section  was  used.  Note  that,  if  desired,  the  additional  weight  of  journal 
sections  may  be  compensated  for  without  increasing  the  number  of 
stations  by  adding  the  weight  as  concentrated  load  at  the  center  of  the 


bearing,  but  such  a load  would  not  be  included  in  the  bearing  reactions 


given  as  output.  All  concentrated  weights,  however,  are  included  in  the 
total  weight  of  the  shaft. 

3.  Where  the  thrust  collar  is  forward  of  the  low  speed  gear  and  not 
integral,  its  weight  is  treated  as  a concentrated  load,  since  it  does  not 
add  significantly  to  the  stiffness  of  the  system.  In  Figure  7,  the  thrust 
bearing  is  aft  of  the  low  speed  gear  and  the  thrust  collar  is  integral  with 
the  shaft;  therefore,  it  is  treated  as  part  of  the  shaft. 

4.  a.  Since  the  gear  is  not  a solid  section  (rim,  webs,  etc. ),  an 
equivalent  outside  diameter  must  be  used  to  obtain  its  proper  lateral 
stiffness.  The  outside  diameter  of  the  hub  section  may  be  used  as  the 
equivalent  outside  diameter. 

b.  The  difference  between  the  total  weight  of  the  actual  low 
speed  gear  and  shaft  assembly  and  the  calculated  weight  of  the  equivalent 
assembly  is  applied  as  a concentrated  load  on  this  section  so  that  the 
total  weight  will  be  correct. 

c.  In  Figure  7,  this  difference  is  applied  as  one  concentrated 
load  at  the  centerline  of  the  gear,  midway  between  fore  and  aft  bearings. 
Present  practice  is  to  divide  the  concentrated  load  by  the  number  of 
webs,  and  to  apply  these  equal  loads  at  the  points  where  the  webs  meet 
the  hub.  In  most  cases,  either  method  may  be  used. 


1 

5 Bearing  reactions  are  assumed  to  act  at  points  midway  between 
the  fore  and  aft  ends  of  each  bearing.  In  the  case  of  long  bearings  (such 
as  a stern  tube  bearing),  it  may  be  assumed  that  the  load  is  concentrated 
one  shaft  diameter  forward  of  the  aft  end  of  the  bearing  For  a more 
exact  detailed  study,  support  points  at  several  positions  along  a bearing 
may  be  assumed. 

6.  If  the  shaft  is  to  be  studied  with  the  ship  waterborne,  a note 
should  be  made  on  the  sketch  (Figure  7)  to  indicate  what  portion  of  the 
shaft  is  to  be  immersed  A station  must  be  initiated  at  the  beginning  of 
the  immersion. 

7.  The  weight  of  the  propeller  is  applied  as  a concentrated  load 
at  its  center  of  gravity,  if  known;  otherwise,  the  weight  is  considered  to 
act  at  the  middle  of  the  tapered  section.  The  dry  weight  of  an  immersed 
propeller  should  be  reduced  to  take  into  account  the  buoyancy  effect  of 
salt  water  by  the  relationship 

Weight  (water)  = Weight  (air)  x (Density  (prop)- Density  (water)) 

Density  (prop) 

For  bronze  with  a density  of  0.  31399  lb/cu  in.  and  salt  water  with  a 


DATA  PREPARATION 


1.  MGE2  Title  Card:  Columns  1-72  alphanumeric  information  to 
identify  the  case. 

2.  One  card  for  each  station  (maximum  of  400)  in  accordance  with 
FORMAT  (13,  5F12.4,  212,  211). 

Column  No  Explanation  

1-3  Station  number  right  justified  (AM,  AA2, ...  , 105,  106, 
etc. ).  Last  station  must  have  9 in  column  1 as  sentinel. 

*4-15  Distance  of  beginning  of  station  from  Station  1,  in  inches. 


*16-27 

*28-39 

*40-51 

*52-63 

64  blank 

65 

66  blank 

67 


Outer  diameters,  in  inches. 

Second  diameter,  in  inches.  Inner  diameters  of  bushing  or 
sleeve  if  used,  or  diameters  or  bore  if  hollow  shaft. 

Third  diameter,  in  inches.  Diameter  of  bore  in  a hollow 
shaft  which  also  has  bushing  or  sleeve. 

Concentrated  weight  (if  present)  minus  values  for  down- 
ward weight. 


♦♦Material  of  outer  diameters 


1,  if  steel 

2,  if  bronze 


♦♦Material  of  inner  diameter 


1,  if  steel 

2,  if  bronze 


68  Bearing  position  - blank,  if  no  support;  1,  if  bearing 
present. 

69  1,  if  shaft  in  water;  blank,  if  in  air 

70  1,  if  sand  bore;  blank,  if  no  sand. 


* The  above  numbers  may  be  punched  anywhere  in  field,  and  the 
decimal  point  must  be  punched. 

**  If  material  is  other  than  steel  or  bronze,  new  constants  for  density 
and  modulus  of  elasticity  must  be  given. 


MEG5  - BEARINGS  RAISED  OR  LOWERED 


PROGRAM  DESCRIPTION 

The  table  of  influence  numbers  obtained  from  MGE2  enables 
the  calculation  of  reactions  when  certain  bearings  are  moved  vertically 
off  the  straight  line.  However,  MGE5  must  be  used  to  find  the  shear 
force,  being  moment,  slope,  and  deflection  at  each  station  caused  by 
raising  or  lowering  certain  bearings 

When  a study  of  shaft  characteristics  from  MGE2  suggests 
that  a better  alignment  could  be  achieved  by  raising  or  lowering  certain 
bearings,  MGE5  should  be  run  to  determine  whether  the  repositioning 
imposes  excessive  shear  forces  or  bending  moments  on  the  shaft. 

The  same  general  computational  procedure  is  used  for  both  MGE5 
and  MGE2.  However,  since  the  influence  numbers  and  weight  and 
stiffness  factors  of  the  shaft  system  are  not  affected  by  changing  the 
elevation  of  bearings,  the  output  of  this  MGE5  program  does  not  list  these 
factors,  It  gives  instead  a single  page  of  shear,  moment,  slope, 
deflection,  and  bearing  reactions. 

Restrictions  on  input  are  the  same  as  those  for  MGE2. 


OPERATING  INSTRUCTIONS 


1,  The  program  run  with  7090  IBSYS  Version  13  includes: 

$ JOB 

$EXECUTE  IBJOB 

$IBJOB  FIOCS 

decks  MGE5 
RED5 

Beam  and  rest  of  decks  same  as  MGE2 

$DATA 

MGE5  heading  card  for  input  case  1 
One  to  give  bearing  change  cards 
MGE2  Reading  card  for  input  case  1 
One  card  per  station  for  input  case  1 
One  card  for  each  constant  change  (omit  if  no  material 
constant  changes) 

Sentinel  card  with  column  2 punched 

Additional  sets  of  input  data  if  called  for  by  sentinel 


Final  sentinel  has  9-or  7-punch  in  column  2. 
i END  OF  FILE 

2«  Output  of  one  page  per  case  is  written  on  SYSOU1  for  program 
control  printing.  Plot  option  is  as  in  MGE2. 

3.  Execution  time  is  similar  to  MGE2. 

4.  The  program  run  on  LARC  includes: 

Data  cards  for  Program  MGE5,  followed  by  six  cards  with 
END  OF  TAPE  punched  in  columns  1-12,  are  converted  to  tape  for 
mounting  on  Unit  11. 
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Object  code  tape  containing  instructions  for  routines 


MGE5 

READ5 

BEAM 

DEFL 

INOU 

INFN 

AM  MAT2L 


is  mounted  on  Unit  13. 


Output  is  written  on  Tape.  20  for  1/1  LSC  loop  printing. 


; * 4/’> 


DATA  PREPARATION 


1.  MGE5  Title  Card:  Columns  1-72  alphamaric  information  to 
identify  how  case  differs  from  similar  MGE2  run. 

2.  One  to  five  bearing  change  cards  using  FORMAT  (II,  5F11.  8). 


Column  No.  Explanation 


13-23 

24-34 

35-45 

46-56 


Blank  or  0,  if  not  last  card:  9,  if  last  change  card. 
Relocation  of  bearing  1 from  straight  line. 
Relocation  of  bearing  2 from  straight  line. 
Relocation  of  bearing  3 from  straight  line. 
Relocation  of  bearing  4 from  straight  line. 
Relocation  of  bearing  5 from  straight,  line. 


Next  card  would  contain  relocation  for  next,  five  bearings.  Values 
are  given  in  inches;  location  of  bearing  may  be  punched  anywhere  in 
field  and  the  decimal  point  must  be  punched. 


Plus  value  is  bearing  above  the  line;  minus  value  is  below. 

3.  A complete  set  of  MGE2  input  cards,  including  case  sentinel. 
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APPENDIX  C 

NOTES  ON  MODIFICATION  FOR  OTHER  COMPUTER  SYSTEMS 

1.  A programmed  double-precision  arithmetic  was  used 
internally  in  matrix  inversion  on  7090.  A computer  with  greater 
significance  might  use  a single-precision  version,  or  the  LARC 
version  MAT2L  might  be  modified  to  avoid  the  assembly  language 
subroutine. 

2.  Use  of  SC  4020  option  assumes  a library  of  MAP  coded 
routines  available  on  NSRDC  system  tape-  A FAP  version  of  these 
routines  exists 

3 For  FORTRAN  II  use  insert  F on  library  and  open  subroutine 
names  as  follows: 

MGE2  0104,  0108,  0109 

MAT2  0040,  0037 

A FAP  version  of  DPAF  also  exists  Card  MGE2  0159  may  be  changed 
to  reflect  standard  system  stop;  e.  g.  , CALL  EXIT. 

4=  Total  allowed  number  of  stations  (400)  would  handle  an 
aircraft  carrier.  For  a smaller  computer,  all  dimensions  of  400,  401 
802,  and  1620  may  be  reduced  equivalently. 

Card  RED2  0032 
or  RED5  0032 

must  reflect  the  maximum  allowed  number  of  stations. 
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